Nonunitary geometric phases: a qubit coupled to an environment with random noise 
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We describe the decoherence process induced on a two-level quantum system in direct interaction 
with a non-equilibrium environment. The non-equilibrium feature is represented by a non-stationary 
random function corresponding to the fluctuating transition frequency between two quantum states 
coupled to the surroundings. In this framework, we compute the decoherence factors which have 
a characteristic "dip" related to the initial phases of the bath modes. We therefore study differ- 
ent types of environments, namely ohmic and supra-ohmic. These environments present different 
decoherence time-scales than the thermal environment we used to study. As a consequence, we com- 
pute analytically and numerically the non-unitary geometric phase for the qubit in a quasi-cyclic 
evolution under the presence of these particular non-equilibrium environments. We show in which 
cases decoherence effects could, in principle, be controlled in order to perform a measurement of the 
geometric phase using standard procedures. 

PACS numbers: 05.40.-a;05.40.Ca;03.65.Yz 



I. INTRODUCTION 

The spin-boson model is studied in a variety of fields, 
such as condensed matter physics, quantum optics, phys- 
ical chemistry and quantum information science [l[ in or- 
der to describe non-unitary effects induced in quantum 
systems due to a coupling with an external environment. 
For a quantum system, the influence of the surroundings 
plays a role at a fundamental level. When the environ- 
ment is taken into consideration, the system dynamics 
can no longer be described in terms of pure quantum 
states and unitary evolution. From a practical point of 
view, all real systems interact with an environment to 
a greater or lesser extent, which means that we expect 
their quantum evolution to be plagued by non-unitary 
effects, namely dissipation and decoherence. 

Most theoretical investigations of how the system is af- 
fected by the presence of an environment have been done 
using a thermal reservoir, usually assuming Markovian 
statistical properties and defined bath correlations 0, Q 
(there are also works on non-Markovian models as, just 
for example, Q). However, there has been some growing 
interest in modelling more realistic environments, some- 
times called "composite" environments [5 7]. In fact, 
the are many situations where the environment is bet- 
ter modelled by a non-equilibrium bath. Quantum dy- 
namics in non-equilibrium environments has been previ- 
ously considered by some recent investigations. For ex- 
ample, light-induced ultra-fast coherent electronic pro- 
cesses in chemical or biological systems may occur on 
sufficiently short time scales Q. In these cases, initial 
non-equilibrium states induced in the bath through the 
coupling among system and environment, might not have 
the chance to reach equilibrium rapidly. Then, the tran- 
sient non-equilibrium bath dynamics may undergo a non- 
trivial interaction with the system of interest in compa- 
rable time scales. Gordon et al discussed the control of 
quantum coherence and the suppression of dephasing by 



stochastic control fields @. In Q, the decoherence pro- 
cess induced by a non-equilibrium environment described 
by several equilibrium baths at different temperatures, is 
discussed. Therein, it was suggested that the effect of 
such environment on the quantum system could be de- 
scribed as the effect done by a single effective bath with a 
time-dependent temperature. The decoherence of single 
trapped ions coupled to engineered reservoirs, where the 
internal state and coupling can be controlled was studied 
in 0. 

In this context, we shall describe a simple model which 
gives a different insight into the behaviour of a quantum 
system coupled to an environment that is not at thermal 
equilibrium. Herein, we study the dynamics of quantum 
coherence in non-equilibrium. We consider a two level 
quantum system in a non-equilibrium bath, represented 
by random perturbations with non-stationary statistics. 
Therefore, we shall study how the quantum system is 
affected by the decoherence induced by the environment. 
We shall compare this decoherence process with the usual 
results for a thermal environment. 

From another point of view, a system can retain the 
information of its motion when it undergoes a cyclic evo- 
lution in the form of a geometric phase (GP), which was 
first put forward by Pancharatman in optics and later 
studied explicitly by Berry in a general quantal system 
[lU . Since then, great progress has been achieved in this 
field. As an important evolvement, the application of 
the geometric phase has been proposed in many fields, 
such as the geometric quantum computation. Due to its 
global properties, the geometric phase is propitious to 
construct fault tolerant quantum gates. In this line of 
work, many physical systems have been investigated to 
realize geometric quantum computation, such as NMR 
(Nuclear Magnetic Resonance) [13[ , Josephson junction 
, Ion trap [lj| and semiconductor quantum dots . 
The quantum computation scheme for the GP has been 
proposed based on the Abelian or non-Abclian geomet- 
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ric concepts, and the GP has been shown to be robust 
against faults in the presence of some kind of external 
noise due to the geometric nature of Berry phase [13- 
[l9j . It was therefore seen that the interactions play an 
important role for the realization of some specific opera- 
tions. As the gates operate slowly compared to the dy- 
namical time scale, they become vulnerable to open sys- 
tem effects and parameters fluctuations that may lead 
to a loss of coherence. Consequently, study of the GP 
was soon extended to open quantum systems. Follow- 
ing this idea, many authors have analysed the correction 
to the GP under the influence of an external thermal 
environment using different approaches (see [2qM23| and 
references therein). 

In this paper, we shall study how the GP is affected by 
the presence of a non-equilibrium environment. We shall 
consider a two-state quantum system coupled to such an 
environment and derive the corresponding decoherence 
factor in SecHU We shall analyse the decoherence pro- 
cess for ohmic and non-ohmic environments. In Sec lIIIl 
we shall derive the GP for a non-unitary evolution of the 
quantum system in the presence of a non-equilibrium en- 
vironment and compute how the GP is corrected in each 
case. Finally, in Sec lIVi we shall make our final remarks. 



II. PURELY DEPHASING SOLVABLE 
SPIN-BOSON MODEL 

A paradigmatic model of open quantum systems is a 
two-state quantum system coupled to a thermal environ- 
ment. This is a particular case of the spin-boson model 
by A. Leggett [TJ (where the tunnelling bare matrix el- 
ement is A = 0) and has been used by many authors 
to model decoherence in quantum computers (2~H and, 
in particular, it is extremely relevant to the proposal for 
observing GPs in a superconducting nano-circuit [25| . In 
spite of its simplicity, this model captures many of the 
elements of decoherence theories and sheds some insight 
into the modification of the GPs due to the presence of 
the environment. The interaction between the two-state 
system and the environment is entirely represented by a 
Hamiltonian in which the coupling is only through cr z . In 
this particular case, [er z ,iJi n t] = and the correspond- 
ing master equation for the reduced density matrix, is 
much simplified, with no frequency renormalization and 
dissipation effects. In other words, the model describes 
a purely decohering (dephasing) mechanism, solely con- 
taining the diffusion term T>(t) [2l[. In such a case, it is 
easy to check that 

p roi (t)=e- nt e- A W ftoi (0), 

is the solution for the off-diagonal terms of the reduced 
density matrix (while populations remain constant) and 
A = J dsV(s). ft refers to the angular frequency of 
precession of a spin precessing the z axis as ruled by the 
isolated from the environment Hamiltonian Ho = \^1<j z 
(responsible for the unitary evolution). The spin-boson 



model is the one used in Ref.[2l| in order to present a 
solvable model to study how the GPs are corrected by de- 
coherence in open systems. In that framework, we have 
studied not only how the GPs are corrected by the pres- 
ence of the different type of environments but estimated 
the corresponding times at which decoherence become ef- 
fective as well. These estimations should be taken into 
account when planning experimental setups, as the one 
performed in Ref. [33|, where using a NMR quantum 
simulator, the geometric phase of an open system un- 
dergoing nonunitary evolution han been obtained. The 
GP was computed in a tomographic manner, measuring 
the off-diagonal elements of the reduced density matrix 
of the system. This study of the GP in the nonunitary 
regime is particularly important for the application of 
fault-tolerant quantum computation (see |34j as an ex- 
ample of measuring the Berry phase in a solid-state qubit 
where there is an important geometric contribution to de- 
phasing that occurs when geometric operations are car- 
ried out in the presence of low- frequency noise). 

In the present paper, we shall adopt a different model 
of decoherence than the one in 211. We are concerned 



with non-equilibrium situations, in which the qubit (the 
main quantum system) is coupled to a non-equilibrium 
bath. The two-level quantum system presents an energy 
gap E2(t) —Ei(t) = hLo(t) which fluctuates due to the in- 
fluence of the environment, where Ej(t), with j = 1,2 is 
the instantaneous energy of state j as perturbed by the 
surroundings. Following the idea proposed in [26j], the 
bath is represented by a random function of time cor- 
responding to the transition frequency of the two-state 
quantum system. In contrast to the usual treatment, the 
statistical properties of this random function arc non- 
stationary, corresponding physically to, for example, im- 
pulsively excited phonons of the environment with initial 
phases that are not random, but which have defined val- 
ues at t = 0. Due to this assumptions, this environment 
is not at thermal equilibrium. The time-dependent fre- 
quency is written in the form Lj{t) = fl + 5u(t), where 
5ui(t) is defined as 



5w{t) = ^2 c k cos(uj k t + 9 k {t)). 



(1) 



fc=i 



The Fourier components c k are positive constants re- 
lated to the spectral density of the environment and the 
coupling of the bath modes to the system. It is impor- 
tant to mention, that, in this model, the randomness en- 
ters through the non-stationary distribution of random 
phases 9 k (t), which are given by 9k{t) = 9 k (0) + x k {t). 
The random function x k (t) satisfies a diffusion equation 



d t P k (x,t) = D k d 2 x P k {x,t) 



(2) 



where P k (x,t) is a time-dependent probability distribu- 
tion and D k is the diffusion constant. The quantity x k 
is an angle, so P(x,t) is a function with period 2ir. The 
time-dependent probability distribution for component 
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k that solves Eq. (|2j) with an initial localized condition 
P(x,0) = S(x) is 



1 



P k ( x ,t) = — + Z 

Z7T ^ 



1 e -n 2 D k t 



,fe 'cos(nx). (3) 



This means that, physically, the phase of each component 
of the random force is not random at t = 0, when an 
impulsive excitation creates a quantum coherence in the 
system, but decays to an uniform l/2n distribution under 
diffusive evolution with diffusion constant Dk (26j . 
Following this approach, it is possible to evaluate, 
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the solution for the off-diagonal element of the density 
matrix (while the populations remain constant again). 
Here, we denote with < ... > the non-equilibrium average 
over the non-stationary random bath and T(t) is defined 
as the decoherence factor. 

By considering the typical factor fk(t) = 
cxp ( — i J Q Sujk(s)ds), one has to do some algebra to 
obtain |.F(t)|. This mainly consists on performing the 
time integral and the ave rag ing fk{t) over the distribu- 
tion probability Pk(xk, t) [26[. After these computations, 
a simple but accurate approximation can be obtained, 
namely 



TOi=in^(*)i^ e " 



(5) 



with 



p{t) = - / du,I(w) 



(6) 



f - e- Ant + [e- zut - e~ wt ) cos(2(wt + 0(w))) 



ft is important to note that in Eq.©, the continuum 
limit has already been taken (in the number of bath 
modes) and the diffusion constant has been assumed 
D(ui) = D for simplicity. 

In order to study the decoherence process induced in 
the system by the presence of a non-equilibrium environ- 
ment, we define a widely usedphysical spectral density 
I(ui) = 4- f /A 2 uj n /A n - 1 e- LU / A Qj, where 7 is the dimen- 
sionless dissipative constant and A is the cutoff frequency. 
On general grounds, A is the biggest frequency present in 
the environment, i.e. the frequency range of the environ- 
mental modes. In particular, the case with n = 1 is the 
"ohmic" case and the one with n > 1 is the "supraohmic" 
case. The ohmic environment is the most studied case 
in the literature, for example in the quantum Brownian 
motion paradigm, and produces a dissipative force that 
in the limit of the frequency cutoff A — > is propor- 
tional to the velocity. The supraohmic case, on the one 
hand, is generally used to model the interaction between 
defects and phonons in metals [I[ and also to mimic the 



interaction between a charge and its own electromagnetic 
field (see for example J27| )■ In particular, the use of the 
supra-ohmic case can be used as a toy model to study 
decoherence process in quantum field theory [28j . 

This model of non-equilibrium is characterized by a key 
quantity which considers the effect of the initial phases 
of the bath modes in the function 0{u>). In this case, we 
consider a linear dependence such as 6{uS) = It is 

interesting to have the possibility to control dephasing 
by varying the single parameter A. Following Eq.©, the 
decoherence factors can be exactly calculated and they 
are given by 
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l-4A 2 (t-A) 2 
~r (l+4A^(t-A)^)a 



(7) 



for the ohmic case, and 
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-67 e - 4£ "( 1 1 c 2Dt )\c 2Dt 1 i- 24 ^-^ A ^+ 16 »-*r A ^ 



(l+4(i-J) ! A 2 ) 4 



(8) 

for the supraohmic case (we use n = 3 all along this 
article) . 

It is interesting to analyse the asymptotic behaviour 
of the function /3(f). Both types of non-equilibrium envi- 
ronments produce a linear time-dependence for the very 
short time-scale Dt, At <C 1, which induces a decoher- 
ence factor of the form J- ~ exp[— wyt] (where a is a 
constant with proper units). This is similar to the de- 
coherence factor calculated in Ref. [2l[ for the case of 
an ohmic finite temperature environment (just assuming 
that 7 ^ JqKbT). In this limit, decoherence is always 
an efficient process, unless the dissipative constant 7 is 
very small. In the long time limit Dt, At ~S> 1, both 
j3 functions (the one for the ohmic, and the correspond- 
ing to the supraohmic cases) acquire a constant value 
(different for each type of environment). In this long 
time regime, decoherence factor in the ohmic case be- 
haves as T ~ exp[— 7], similar to the decoherence factor 
for the equilibrium supraohmic environment at zero tem- 
perature [2l[. Meanwhile, in the supraohmic case, the 
decoherence factor approaches a long time value given 
by T ~ exp[— 67]. Again, as we mentioned before, in the 
case of small 7, decoherence never occurs, even at very 
long times. Intermediate times are ruled by the specific 
randomness introduced into the model. All in all, it is 
important to note the richness of the model which guar- 
antees known and unknown decoherence processes by the 
correct tuning of the parameters. 

In Fig. [I] we present the behaviour of the decoherence 
factor for a strong dissipative case for both environments. 
As expected, the decoherence factor decays from unity 
to an asymptotic value. Therein, we also present the 
behaviour of the decoherence factor found in , where 
a different spectral density to describe the environment 
has been used. The parameters used in the Figure are 
similar to those used in @ in order to do a better com- 
parison and analysis. Unlike typical studies using the 
master equation in the weak coupling limit (7 ~ c^), in 
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FIG. 1. (Color online) Evolution in time of the decoherence 
factor J-(t) for different models of the environment in the 
strong coupling limit. The red short-dashed line is for an 
ohmic non-equilibrium bath and the long-dashed black line is 
for the supraohmic. For Dt, At 3> 1, the decoherence factor 
behaves as J r hmic(i) ~ exp(— 7) and T sup ra(t) ~ exp(-67). 
We also present a solid blue line for the assumptions made in 
H, a non-equilibrium bath with a Gaussian spectral density. 
Parameters used: 7 = 3; A/ft = 1; ftA = 1; D/Q = 0.5; and 
71 = 3 for the supraohmic environment. 



the present approach there is no constraint for the value 
of 7. Therefore we can use either a strong or weak cou- 
pling as a value for 7. In the Figure, we can note three 
different lines: the solid blue line for the results in Q, the 
red short-dashed line for our ohmic environment and the 
solid black long-dashed line for the supra-ohmic environ- 
ment. Then, it is easy to see that the ohmic case is very 
similar to the one obtained in Q , where a Gaussian spec- 
tral density was considered. In both cases, decoherence 
is very efficient, as expected since we are considering the 
overdamped case (7 > 1). There is also an interesting 
fact: the supraohmic decoherence factor has a smaller 
decoherence time-scale than the other two decoherence 
factors herein considered. This is unlike the case for equi- 
librium supraohmic environment, where decoherence is 
effective only at high temperature. This modelling of the 
environment gives a decoherence factor which drops from 
its initial value toward an asymptotic value {T{t — > 00)) 
after the intermediate time t = A. At this time the sys- 
tem re-phases back to the slowly decaying envelope not 
purely exponential. As A becomes large and positive, the 
decay approaches the envelope function without the non- 
monotonic dip (that occurs at t = A). Non-exponential 
behaviour in the decay of quantum coherence has been 
observed in full many-particle simulations of quantum co- 
herent dynamics under non-equilibrium conditions [29j . 

In order to have a better view of the dip where the 
system "recoheres" for a while, we present the behaviour 
of the decoherence factors for the weak coupling case in 
FigEI Therein, 7 has a smaller value, comparable to 
those we used when dealing with environments in ther- 
mal equilibrium in the underdamped case [2ll [30L l3lj . 
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FIG. 2. (Color online) Evolution in time of the decoherence 
factor F(t) for different models of the environment in the weak 
coupling limit. The red short-dashed line is for an ohmic non- 
equilibrium bath and the black solid line is for the supraohmic 
case. We also present a solid blue line for the assumptions 
made in @], a non-equilibrium bath with a Gaussian spectral 
density (different with respect to the normally used in the 
theory of quantum open systems theory). Parameters used: 
7 = 0.5; A/ft = 1; ftA = 1; D/Q = 0.1. 
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FIG. 3. (Color online) A more complex modelling of the initial 
phases of the bath modes by considering 6(uS) = —\u> 2 in 
an ohmic (blue solid line) and supraohmic (black solid line) 
environment. Parameters used: 7 = 3; A/ft = 1; ftA = 1; 
D/Q = 0.1. 



The dip is obtained by introducing, as we mentioned 
before, the simple relation 9{uS] = —\oj for the initial 
phases of the bath modes in the modelling of the envi- 
ronment. Even though this assumption is a minimalistic 
model, it allows to have some kind of control in the de- 
coherence process which in turn can be useful in experi- 
mental setups where decoherence is always an obstacle to 
overcome. This result agrees with the one in Q, where 
it was shown that non-equilibrium decoherence can be 
slowed down in a controlled manner as compared to the 
corresponding equilibrium situation. 

A different modelling of the initial phases of the bath 
modes can, in principle, be adopted. However, herein we 
use the linear one just for simplicity. A complex assump- 
tion can be, for example, 9(lu) = —\oj 2 . The derivation 
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of the decoherence factor is somewhat more difficult and 
is not worth writing explicitly here. Anyway, the deco- 
herence factor for a quadratic behaviour in oj is presented 
in Fig. [3] for an ohmic and supra-ohmic non-equilibrium 
environment. In such a case, it is important to note a 
more complicated structure of dips in the decoherence 
factor. 



III. APPLICATION: GEOMETRIC PHASE FOR 
A QUBIT COUPLED TO A NON-EQUILIBRIUM 
ENVIRONMENT 

In order to compute the GP and note how it is cor- 
rected by the environment, we shall briefly review the 
way the geometric phase can be computed for a system 
under the influence of external conditions such as an ex- 
ternal bath. In Ref. [Iflj , a quantum kinematic approach 
was proposed and the geometric phase (GP) for a mixed 
state under non-unitary evolution has been defined as 

4> G = arg{ ]T vM0)e*(T)<**(0)|tt fc (T)) * 



-J T «"<**! 



(9) 



where Ek(t) are the eigenvalues and the eigenstates 
of the reduced density matrix p r (obtained after tracing 
over the reservoir degrees of freedom). In the last defini- 
tion, r denotes a time after the total system completes a 
cyclic evolution when it is isolated from the environment. 
Taking into account the effect of the environment, the 
system no longer undergoes a cyclic evolution. However, 
we shall consider a quasi cyclic path V : t e [0, t], with 
r = 2tt/CI (ft is the system's characteristic frequency). 
When the system is open, the original GP, i.e. the one 
that would have been obtained if the system had been 
closed <f>u, is modified. This means, in a general case, 
that the phase can be interpreted as cf>Q = (/>u + ^g, 
where S(f>G depends on the kind of environment coupled 
to the main system @, [HM [H [H] . 

Assuming an initial quantum state of the form 

|V(0)>=cos(|)|0>+ sin(|)|l>, (10) 
its evolution at a later time t, is 

|V(t) >= e~ lUt cos((9 + )|0 > + sin(6>+)|l >, (11) 
where 

sm(6 )\F(t)\ 



cos(# + ) 



sin(0+) 



sm 2 (9 )\F(t)\ 2 + 4(£+ _ cos 2(^))2 ' 



2( £+ -cos 2 (^)) 
yJsm 2 (6 )\T(t)\* + 4(e+ - co S 2(f ))* 



(12) 



(13) 



and £ + the eigenvalue of the reduced density matrix, 
namely 



cos 2 (0 o ) + sin 2 (0 o )|^)P 



(14) 



while e_ does not contribute to the geometric case since 
E-(t = 0) = (see definition En.©). 

As in our previous works HH HH, [33j], the GP is 
obtained by computing eigenvectors and eigenvalues of 
the reduced density matrix and using Eq.©, 



n 



dt. 



(15) 



In Figs. 0] and [S]we plot the environmentally induced 
correction to the unitary phase, \84>g\ (normalized by 
the value of 4>u) as a function of the system's initial 
quantum state (6q) and the dissipation induced in the 
quantum subsystem due to the presence of the random 
environment (7). In both Figures we have considered a 




FIG. 4. (Color online) Behaviour of the geometric phase 
as a function of the initial state of the quantum system do 
(measured in radians) and the dissipation of the environment 
(dimensionless 7) induced by an ohmic non-equilibrium en- 
vironment in a cycle. Parameters used: A/Q = 1; S1A = 1; 
D/Q = 0.1. 

wide range of values for 7, considering both weak and 
strong coupling between system and environment. For 
small values of 7, the GP behaves very similarly to the 
unitary GP which is <pu = 7r(l + cos#o)- However, as we 
increase the value of 7, there is a notable change in the 
curvature as a function of 9o , leading to more values of 9q 
with a null GP. As expected, the more decohering envi- 
ronment, the less survival of the GP. In agreement with 
Fig 12 we can see that the geometric phase is less cor- 
rected (with respect to the isolated case) in the presence 
of an ohmic non-equilibrium environment rather than of 
a supraohmic one. This can be noted by the fact that 
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FIG. 5. (Color online) Behaviour of the GP as a function of 
the initial state (do in radians) of the quantum system and 
the dissipation of the environment (dimensionless 7) induced 
by a supraohmic non-equilibrium environment in a cycle. Pa- 
rameters used: A/O = 1; f2A = 1; D/Q = 0.1. 

Fig |4] remains a smooth function of #0 for bigger values 
of 7 than FigJSJ in which case the phase rapidly behaves 
different as a function of the dissipation constant. 

The GP cannot be fully computed analytically but we 
can perform an expansion in powers of the coupling con- 
stant, to obtain an accurate approximation of it [2l|-|23J. 
Hence, we expand in powers of 7 the cos 2 9 + in Eq. (|15l) . 
using the definition of the decoherence factors for each 
environment, namely Eqs.Q and (jHJ. As mentioned be- 
fore, the correction to the GP is defined as 54>c while 
(j>u is the unitary GP. In the case of the ohmic non- 
equilibrium environment, the correction to the unitary 
GP is given by, 

Hg„ =1 ~ 7r7sin 2 (6>o)cos(6> ) 

+ 7^e- 2m sin 2 (# o )cos(0 o ). (16) 

Proceeding the same way for the supraohmic environ- 
ment, the correction to the unitary GP is 

S(j)c n=3 ~ 671-7 sin 2 (#o) cos(#o) 

+ 7^ e - 2m sin 2 (tfo)cos(fl ). (17) 

The corrections of the GP for both environments agree 
for small values of 7 in Figs. @] and [5j In both cases, 
the dependence with the parameter A is exponentially 
negligible. Another interesting feature of the corrections 
of the GP is that they depend on the initial angle of 
the quantum state, and this dependence is in agreement 
with the ones obtained for a two-level system in interac- 
tion with environments at equilibrium 0, I2H - I23I [H, HH . 
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FIG. 6. (Color online) Comparison between exact GP (red 
dashed line), in the presence of an ohmic environment, and 
the first order perturbative correction (blue solid line) from 
Eq. (|16|) . Perturbative result is in good agreement with the 
exact result for a long range of values of 7. Parameters used: 
D/Q, = 1, OA = 1, A/O = 1. 

Neglecting the small correction induced by A (which is a 
correct assumption seen Figs. S]and[5]), both cases are 
similar to the very low temperature corrections found in 
[2l| for the case of thermal environments. In Figs. [6] and 
we show the range of validity of the first order per- 
turbative expansion in powers of 7. In FigEl it is clear 
that the perturbative result (solid line) of Eq. ([T6"f is in 
excellent agreement with the exact result (dashed line), 
even for not too small values of the coupling strength 
parameter 7. Fig. [7] shows that Eq . ([]"?]) is also a good 
approximation to the exact result (dashed line), but only 
for very small values of 7. 




FIG. 7. (Color online) Comparison between exact geometric 
phase (red dashed line) , in the presence of an supraohmic en- 
vironment, and the first order perturbative correction (blue 
solid line) from Eq. (|17[l . Perturbative result is a good ap- 
proximation for really small values of 7. Parameters used: 
D/Q, = 1, OA = 1, A/O = 1 

It is an interesting feature to study which is the in- 
fluence of the observed dip (in the decoherence factor) in 
the behaviour of the GP. In this model of non-equilibrium 
environment, the parameter A, which enters through the 
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FIG. 8. (Color online) Behaviour of the geometric phase as 
a function of the initial state 9o (in radians) of the quantum 
system and the initial phases of the bath modes (A) in a cycle 
for an ohmic non-equilibrium environment. Parameters used: 
A/Q = 1; 7 = 3; D/Q = 0.1 

random phases, sets the position at which the "dip" or 
"recoherence" takes place. As we have seen, there is 
a time-scale when the system seems to gain coherence 
(i ~ A). In FigEl we present the correction to the GP 
\8<pG I as a function of the initial state of the quantum sys- 
tem (6*o ) and the initial phases of the bath modes (OA)for 
an ohmic non-equilibrium environment. Therein, we can 
observe that the geometric phase is not affected by the 
dip in the decoherence factor, as it is computed over a 
quasi-cyclic evolution. In Figj8]we can note that the GP 
has a monotonous behaviour as a function of A. The an- 
alytical estimation of the influence of A in the correction 
to the GP, made in Eqs. (fT!)|) and p7|) . is also checked 
in Fig. El 

IV. FINAL REMARKS 

The geometric phase of quantum states is an issue 
worth of attention. It could be a potential application 
in holonomic quantum computation since the study of 
spin systems effectively allows us to contemplate the de- 
sign of a solid state quantum computer. However, deco- 
herence is the main obstacle to overcome. All realistic 
quantum systems are coupled to their surroundings to 
a greater or lesser extent. Furthermore, in most cases 
of practical interest, quantum systems are subjected to 
many noise sources with different amplitudes and corre- 
lation times, corresponding de facto to a non-equilibrium 
environment. 

Herein, we have presented a simple case to illus- 
trate the general phenomenon of dephasing in a non- 



equilibrium bath. We have studied the decoherence pro- 
cess of a quantum system in interaction with an initially 
non-equilibrium bath that can be controlled by manipu- 
lating the nature of the relative initial phases of the bath 
modes. The decoherence factors computed here suggest 
that by engineering these initial phases, the dephasing 
of the subsequent quantum evolution can potentially be 
controlled. Wc have found similarities and differences of 
the decoherence process between the environment pre- 
sented here and thermal environments studied in previ- 
ous works. The model presented here is another proposal 
for engineering reservoirs in a manner reminiscent of a 
coherent control experiment using shaped pulses [HI . In 
this model, the control parameter A is derived not from 
a laser pulse, but rather from well-defined phase rela- 
tions between the modes of the bath. Another possible 
candidate for realizing this decoherence environment, is 
to use the artificially generated fluctuating environments 
with NMR. It could be possible, in principle, to use the 
quantum simulator of Ref [33j to generate the fluctuating 
phase 0(lu) of the present proposal. 

The analysis of the effect produced by decoherence on 
the GP is crucial at the time to design an experimental 
setup to measure the GP using, for example, interfcr- 
ometry. We found that the convenient non-equilibrium 
environment to observe GPs is the weak coupled ohmic 
case. It is important that there is no restriction about 
zero temperature environments in this case, as it was 
found in plj . as the most convenient scenario. In this 
framework, these kind of environments could become 
a proper experimental setup for the observation of 
the geometric phase. Therefore, we have computed 
the geometric phase for an ohmic and supra-ohmic 
non-equilibrium environment and seen how they deviate 
from the unitary geometric phase. So far, we have seen 
that the characteristic dip of the decoherence factor does 
not affect much the geometric phase and that ohmic 
non-equilibrium environment arc not as destructive as 
the supraohmic non-equilibrium ones or the thermal 
environments we are so used to see in the literature. 
The effect done on the geometric phase by the ohmic 
(or/and supraohmic) non-equilibrium environment can 
be seen as similar to the one done by a single reservoir 
with an effective temperature, as the non-equilibrium 
environment model proposed in p|. In the very weak 
coupling limit, we have evaluated the corrections 
induced by the non-equilibrium environment on the 
unitary GP, showing that there is a small (exponentially 
suppressed) correction due to the random parameter A. 
More general models should be analysed in a future work. 
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